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Abstract 

Many studies on microscopic systems deal with Brownian particles embedded in media whose 
densities are different from that of the particles, causing them either to sink or float. The proximity 
to a wall modifies the friction force the particle experiences, being, instead of the Stokes force 
assumed by Einstein and Langevin, a function of the separation distance, as described by Faxen. 
In this work we present a thorough analysis of bidimensional Brownian motion using monodisperse 
suspensions of latex beads with controlled temperatures. By means of Langevin's model and 
Faxen's correction, we calculate the mean separation distance between the spheres and the bottom 
of the confining vessel. 
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I. INTRODUCTION 



In 1827, Brown observed organic and inorganic microparticles having an uninterrupted 
and haphazard motion in fluid media.- With the same result for many materials, he con- 
cluded that all matter should be constituted by a "primitive molecule". For 50 years the 
new phenomenon, later known as Brownian motion, did not received any contribution to 
its explanation in physical nor mathematical terms. In 1877 Wiener demonstrated that the 
motion could not be due to convection currents and Dancer discarded electrostatic effects as 
the reason.- It was Desaulx and Carbonnelle who suggested that Brownian motion should 
constitute a physical phenomenon per se, where molecules of the medium hit the micropar- 
ticle causing it to move randomly. And in 1888 Gouy demonstrated it qualitatively, by 
showing experimentally that Brownian motion is independent of the nature of both, the 
medium and the particle, although it is certainly affected by temperature.- 

In his search for evidence of the same idea, the existence of molecules, Einstein argued in 
1905 that the random motion of microscopic particles must be due to the multiple collisions 
with the yet much smaller molecules of the medium.- He derived his celebrated formula 
relating the macroscopic diffusion coefficient with the temperature and the viscosity of the 
fluid, and with Avogadro number A^^- Moreover, he went further by proposing a simple 
experiment to measure the mean square displacement (MSD) of a micron-sized particle and 
to obtain, in turn, the diffusion coefficient. Some years later, in 1908, Langevin proposed a 
version of Newton equation for a Brownian particle in a fluid from which Einstein's result 
naturally followed. The developments by Einstein and Langevin are now milestones of the 
general theory of Markovian processes.- 

In 1906, inspired by Einstein's proposal, Jean Perrin calculated A^^ according to Ein- 
stein's expression for the MSD. In his experiments, the trajectories of hundreds of gamboge 
and mastic microparticles with radii in the range of 0.21-5.5/im were followed by hand using 
a camera lucidaA^ The particles were embedded in media with different viscosities, being the 
largest value up to 125 times larger than the smallest one. From this statistics he obtained 
a value of A"^ = 71.5 x 10^^ particles per mol (see page 90 of Ref.-), which was consistent 
with the values he obtained based on different studies, and hence this constituted the first 
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experimental demonstration of Einstein's prediction.- Since that moment, Brownian motion 
acquired an important presence in many topics of physics, moving us to explore it on its 
fundamentals. 

Repeating Perrin's classic experiment with modern means has become a very attractive 
endeavor,-""- not only because of its historical importance but also because it is still a ref- 
erence for understanding macroscopic phenomena where the role of the molecular collisions 
is relevant. It is very important to mention, however, that most studies on Brownian mo- 
tion, including Perrin's, involve particles with densities different from that of the medium, 
and therefore, they either float or sink. This fact is not considered by Einstein's nor by 
Langevin's theories. In their models, the only forces are the dissipative one, considered as 
the Stokes drag for a spherical particle, and the stochastic force due to the molecules of 
the medium. As we point out and study in the present article, once the particles sink to 
the bottom of the container, the interaction with the wall gives rise to a correction to the 
viscosity of the medium that substantially modifies Einstein relation. This correction, well 
known by the scientific community working on optical micromanipulation, may be taken 
into account by invoking to Faxen's modeL-^i^ 

In this work the ultimate goal is not to find the value of Avogadro number A^^, but 
to reconsider the approximations assumed to describe typical microparticle systems. Our 
interest is to study Brownian motion taking into account the presence of a boundary. Our 
scheme is the same as Perrin's, but instead of the camera lucida we use a coupled-charged 
device (CCD) and perform a thorough computer analysis of the data. We record the mo- 
tion of individual latex microspheres and from the videos we obtain their positions as a 
function of time, which allows us to calculate the MSD. In order to avoid hydrodynamic 
and electrostatic interactions among particles, we use very low concentrations of particles 
within an ionic solution. The paper is organized as follows: Section [TTl briefly describes the 
Brownian motion model in a timescale larger than the relaxation time of the system, the 
diffusive regime. In Section IIIH the experimental setup and the data analysis algorithm 
are described. At the end, A^i is calculated without the effective viscosity corrections by 
analyzing the MSD as function of time. The comparison of this result with the known 
standard value, shows the necessity of including a different description of the system. In 
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Section ITV] we introduce Faxen's correction and, assuming the standard value of A'^, the 
mean vertical distance of the particles from the bottom is calculated. On the basis of our 
results, we conclude with some remarks and suggestions for further experiments. 

II. BROWNIAN MOTION MATHEMATICAL MODEL 

In Langevin's approach, the motion of a spherical microparticle in thermal equilibrium 
with its surroundings can be expressed by Newton second law where the total force is the 
addition of a dissipative term and a time- dependent stochastic force. ^ The friction force 
is the mechanism by means of which the particle looses heat, while the stochastic force, 
i{t), representing the collisions of the molecules of the medium with the particle, returns 
energy to the particle. One should remark that energy can be transferred just by collisions 
and rotations since the kinetic energy of the molecules is not enough to produce internal 
excitations in the particle. As the translational contribution to the kinetic energy is much 
larger than the rotational contribution,AJi we neglect the rotational dynamics in the model. 

According to Stokes, for a spherical particle falling in an unbounded medium with small 
Reynolds number, the dissipative force is a drag proportional to the velocity of the particle, 
with coefficient 7 = Snrjd, being rj the viscosity of the medium and d the diameter of 
the particle. Meanwhile f(t) is a white noise force whose first moment is null and its 
second moment is f(t)f(t') = a6{t — t'), which means that it is completely uncorrelated 
with itself at different times (this is true for time intervals larger than the correlation time 
of the immersion medium, inherent to the duration of the collisions between molecules). 
The coefficient a = 2RT'j/Na as given by the Fluctuation-Dissipation Theorem, expresses 
the fact that the particle is in thermal equilibrium with the medium; T is the absolute 
temperature and R the gas constant. The components x, y and z of f(t) are completely 
uncorrelated with each other.— Langevin motion equation then reads: 

m^ = -7v + f(t), (1) 

where m is the mass of the microparticle and v its velocity. When the density of the particle, 
Pp, is different from the density of the medium, p^, one should include a third term given 
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by: 



(2) 



with g the gravitational acceleration. F causes in general the particle to remain on the 
bottom or at the top. When the particle is close to a wall, Faxen found that the consequence 
is an increase of the effective viscosity the particle experiences. The drag force, then, is better 
represented by —k'jv, where k is given by:— "i^*^ 
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being h the distance from the center of the particle to the wall. Notice that k tends to 1 as 
h tends to infinity, recovering the Stokes force. 



Because of the fluctuations along z-axis, h is a function of time, nevertheless we can 
analyze just the bidimensional motion on the plane assuming h as an average value of the 
vertical position of the particle. In this case the gravitational force needs no longer to be 
explicit in the motion equation. 



Therefore, we write the Langevin-Faxen bi- dimensional Brownian motion equation for a 
particle close to one wall as: 

dv I 

m^ = -fc7v± + f±(t). (4) 

The factor m/ (k'-f) = tr is known as the relaxation time. For At ^ tr the probability of 
finding the particle with velocity vi at ti = to + is completely independent of vq at to? so 
the motion is in a non-correlated regime. Moreover, once the particle relaxes to equilibrium, 
it has a velocity distribution that must be the Maxwell velocity distribution.— For typical 
Brownian systems is ~ lO^^s, which is much smaller than the time resolution for usual 
imaging systems. For instance, our CCD takes 30fps, hence, we are able to observe only a 
random walker. 



In this regime, the term on the left hand side of Eq. (jl]) is negligible when divided by /c7, 
then we have for the bidimensional position, r: 
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dt k'-) 

Equation ([5]) can be solved as an integral equation, 
solution and considering that r = tq, one finds: 



Taking the time average to the 



r - rol = 0, 

a . . RT 
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ro)2| = 2—-\t = 4— —At, (6) 



because of the statistical properties of f(t). In the limit A; — )■ 1, Eq. ([6]) corresponds to 
Einstein's solution: 



|r2 - r^l = ADt, (7) 
with the diffusion coefficient identified as: 



which is the celebrated Einstein relation.- 



III. EXPERIMENTS 

Figure [1] depicts the experimental setup for studying Brownian microparticles, which was 
basically a standard video microscopy system. It consisted of a 100 x immersion objective 
(N.A. 1.25) and two lenses (/Li=50mm and /L2=100mm) playing the role of the eyepiece. 
The total magnification was approximately 140 x. As a condenser we used a Koehler illu- 
mination system, which provides optimum contrast and resolution, and uniform brightness 
in the illumination field.— The video camera was a CCD with a recording speed of 30fps. 
We restrict our analysis to the motion on the horizontal plane, xy. 



Our samples consisted of latex spheres, with either di = 1.0/im or d2 = 2.0/im in diameter 
(size dispersion 4%)^, suspended in a diluted saline solution. This ionic medium avoids 
electrostatic interactions among particles.— The concentrations were low enough to allow 
the imaging of just one particle per video, avoiding also the hydrodynamic interactions in 
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FIG. 1. Setup for viewing microparticles with Brownian motion. 



this way.-i^ The density of the particles, pp =1.05g/cm^, is higher than that of the hosting 
medium, pm =1.00g/cm^, and thus they settle on the bottom of the sample cell. The sam- 
ple cells were shallow cylindrical chambers with 100/xm depth and 0.5cm in diameter made 
with a dark sticker placed between a microscope slide and a cover slip, both previously 
cleaned with a free rinsing surfactant (LiquiNox) and with acetone, then sealed with epoxy 
resin. The samples were placed on a plate with a temperature control system, and the plate 
was attached to an xyz translational stage, allowing us to choose the particle of interest. 
It is worth to mention that it was imperative to guarantee the horizontal position of the 
sample cell, in order to avoid motion currents on the xy plane due to the gravitational force. 

We recorded Q = 10 videos for each particle size for over 1 minute, which corresponds to 
L = 1800 frames per video. The particles were tracked with a Matlab code which analyzes 
frame by frame finding the centroid of the particle image and recording its two-dimensional 
coordinates as a function of time. The input parameters for this processing are: 1) the 
diameter of the object in pixels, 2) the tolerated image noise associated with the change of 
the image resolution when the particle is out of focus (due to vertical Brownian motion), 
and 3) the maximum displacement from one frame to the next one, in order to avoid the 
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tracked particle to be confused with a different one in case many particles are present. Since 
the distances are expressed in pixels, it is necessary to find the equivalence between pixels 
and micrometers. This was made by measuring in pixels the diameter of the spheres and 
comparing with the size reported by the manufacturer. From this process Q data files were 
generated, these include the two-dimensional positions (fq'^^ r['^\ r2'^'', . . . ,r^^^) separated at 
regular time intervals At = (l/30)s. Here g = 1, 2, . . . , Q denotes the different videos. 



Since we are assuming a diffusive Markovian stationary process, each step of the particle 
should be independent of its initial position. We can take advantage of this fact to perform 
averages with very good statistics by defining relative displacements as a|'''^'' = r'f^^ — r[^^ , 
with / = 0, (L — 1 — Smax) and s = 1, 2, ... , Smax- Therefore, for every time t = sAt we 
calculate the MSD as 
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In our case, we restricted our analysis to the time interval t G [0, 10)s, which implies that 
Smax = 300. Therefore, a total of 15 000 data were used for the calculation of each MSD 
point, 1500 per video, thus making negligible the statistical error.— 

Figure H] shows the plots obtained for the MSD as a function of time for both sizes 
of particles at different temperatures, 20, 40 and 60°C. The linear regression through the 
origin for these functions provided us a value for the diffusion coefficient D according to 
Eqs. ([7]) and (|8]), from which we could obtain Avogadro number A^^. For instance, for di 
and T = 20°C, we obtained Na = 9.584 x lO^^part/mol. Uncertainty on this value can 
be due just to three factors: 1) reported dispersion size of the spheres 2) measurement in 
pixels of the image of the particle 3) accuracy of the temperature controller. In any case, 
they are not enough to explain the large difference with reported values of Na by different 
techniques.— This incite us to use Eq. (|6]) instead of Eq. ([7]). 
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FIG. 2. (Color online) MSD vs time for latex particles with di = l.O^m (a) and d2 = 2.0/im (b) in 
diameter. Dot-dashed lines correspond to T = 20° C, dashed lines to 40° C and solid lines to 60° C. 

IV. FAXEN'S BOUNDARY CONDITIONS 

By analyzing Fig. [2l we found the linear dependence of MSD with time, nevertheless the 
proportionality constant is substantially smaller than the Einstein coefficient, 4D. On the 
other hand, the use of Eq. (jl]) makes possible either to determinate A^^ more accurately by 
knowing the average vertical position of the particle, or to determinate the average vertical 
position of the particle, h, by assuming the standard A^^ value. We explored the second 
option. 

Figure [3] shows the correction factor for the viscosity k as function of h, given by Eq. (|3]), 
for particles of diameters di (dot-dashed line) and d2 (solid line). From the values we found 
for the slopes in Fig. |2] and using Eq. ([6]), we were able to determine the values of k for 
each of the studied particle sizes and temperatures. The dotted line and the dashed line in 
Fig. [3] correspond, respectively, to the values of k found for di and d2 a.t T = 20°C, and the 
associated values of h, are determined graphically from the curves for k vs h. 

The latex particles sank to the bottom and Fig. [3] is consistent with this behavior, since 
it implies that the mean vertical position above the bottom surface is around 0.73/im for 
the smallest particle {di) and 1.52/xm for the largest particle (^2)- 
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FIG. 3. (Color online) k vs h for di (dot-dashed line) and d2 (solid line). Straight lines show the 
mean vertical position of di particles (dotted line) and d2 particles (dashed line), at T = 20° C. 

V. SUMMARY 

With a very simple scheme we have observed Brownian motion for 1 and 2/im-diameter 
latex spheres, avoiding electrostatic and hydrodynamical interactions among them, using 
a ionic medium and a very low concentration. The linear behavior of the mean square 
displacement as function of time was shown for the diffusive regime, in which the sampling 
time is larger than the relaxation time of the system. 

Calculations of A^^ obtained from a thorough statistical analysis and according to 
Langevin and Einstein, are very different from standard values. The reason is that the 
assumption of an unbounded system was not fulfilled in this case, where the particles sink. 
This is the most typical situation in systems involving motion of microparticles. The con- 
sequence of the sedimentation is that the drag force cannot be assumed to be the Stokes 
force, as Langevin did. 

We rather employed the Faxen's correction to the drag, which is a function of the ratio 
between the diameter of the spheres and the average separation distance from the wall. 
This provides a way to measure Na_ in a better approximation or, alternatively, to calculate 
the average vertical position of the Brownian particle if a standard value of Na is assumed, 

10 



either of both with the same simple setup. 



Studying Brownian motion can be made in a relatively simple way in an undergraduate 
laboratory, but it is important to keep in mind that its theory has many subtle issues that 
could lead to in an improper interpretation of results if they are not considered. When 
analyzing Brownian systems one must never forget to guarantee the sample to be in thermal 
equilibrium, to avoid any interaction among particles and to consider boundary conditions 
of the system, where gravity may play an important role. 
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